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We explain how a four-body amplitude analysis of the D decay products in the mode — > (K+K~7r+7r~)D 
is sensitive to the unitarity triangle angle 7. We present results from simulation studies which show that a precision 
on 7 of 15° is achievable with 1000 events and assuming a value of 0.10 for the parameter vb- 



1. Introduction 

A precise measurement of the unitarity trian- 
gle angle 7 is one of the most important goals 
of CP violation experiments. 7 is defined as 
a.rg{-V*^Vud/V*i,Vcd), where Vij are the elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) mix- 
ing matrix. In the Wolfenstein convention [T] 
7=arg{T/*J. 

A class of promising methods to measure 7 ex- 
ists which exploits the interference between the 
amplitudes leading to the decays B^ D"K^ 
and B" D"K- (Figure [T]), where the D° 
and D'' are reconstructed in a common final 
state. This final state may be, for example, a CP 
eigenstate such as K+K~ ('GLW method') [5], 
or a non-CP eigenstate such as K+7r~, which 
can be reached both through a doubly Cabibbo- 
suppressed D° decay and a Cabibbo-favoured I)'^ 
decay ('ADS method') |3j. Recent attention has 
focused on self-conjugate three-body final states, 
in particular D — > Kg7r+7r~ 0. Here a Dalitz anal- 
ysis of the resonant substructure in the Kstt+tt" 
system allows 7 to be extracted [4] . The B-factory 
experiments have used this method to obtain the 
first interesting direct constraints on 7 |6l7j . 

Here we explore the potential of determin- 
ing 7 through a four-body amplitude analysis 
of the D decay products in the mode B* 
(K+K^7r+7r^)D K^. CP studies involving am- 
plitude analyses of four-body systems have been 

'^Here and subsequently D signifies either a D*^ or a C. 



proposed elsewhere [4] , and strategies already ex- 
ist for B^ DK^ approaches exploiting singly- 
Cabibbo suppressed decays |5j . Our method ben- 
efits from a final state that involves only charged 
particles, which makes it particularly suitable for 
experiments at hadron colliders, most notably 
LHCb. 

This paper is organised as follows. In Sec- 
tion [2] we summarise the essential features of 
B^ DK^ decays, and state the present knowl- 
edge of the parameters involved, and of the decay 
D -> K+K-TT+TT". In Section [3] a full model of 
B^ (K+K^7r+7r^)D K^ decays is developed, 
which is then used within a simulation study to 
estimate the precision on 7 which may be ob- 
tained through a four-body amplitude analysis. 
We conclude in Section [H 

2. B± DK± decays and D -> K+K-tt+tt" 

Let us define the amplitudes of the two dia- 
grams illustrated in Figure [T] as follows 

A(B-^D"K-) = Ab, (1) 
A{B- ^B"K-) = Asrse^f^^-^). (2) 

Here the strong phase of As is set to zero by con- 
vention, and Sb is the difference of strong phases 
between the two amplitudes. 7 represents the 
weak phase difference between the amplitudes, 
where contributions to the CKM elements of or- 
der and higher (with A being the sine of the 
Cabibbo angle) have been neglected. In the CP- 
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Figure 1. The diagrams for B~ D°K~ and 
B~ D"K^. There is a relative phase oi 6b — 7 
between the two amphtudes, and a relative mag- 
nitude of rs. 



conjugate transitions 7 — s- —7, whereas 6b re- 
mains unchanged, is the relative magnitude 
of the colour-suppressed D°K~ process to 

the colour-favoured B^ D°K~ transition. Pre- 
liminary indications as to the values of 7, 6b and 
tb come from the B^ DK^,D — > Kstt+tt^ 
analyses performed at the B-factories |6I7| . Fits 
to the ensemble of hadronic flavour data also pro- 
vide indirect constraints on the value of 7 |8l9j . 
These results lead us to assume values of 7 = 60° 
and 6b — 130° for the illustrative sensitivity stud- 
ies presented in Section [3l We set rs to 0.10, 
which is the approximate average of the Dalitz 
results and the lower values favoured by the ADS 
and GLW analyses [TOiTT] . 

Results have recently been reported from an 
amplitude analysis of the decay D K+K^tt+tt^ 
[12j . which shows that the dominant contribu- 
tions come from D AP and D VV modes. 
Earlier measurements of D ^ K+K^tt+tt^ were 



published in [13]. Our sensitivity studies for the 
7 extraction, presented in Section [3l are based on 
the results found in [12]. 

The branching ratio of the mode B^ — > 
(K+K~7r+7r~)D can be estimated as the prod- 
uct of the two meson decays, and found to be 
9.2 X 10"'' [l4j. This channel is particularly well 
matched to the LHCb experiment, on account 
of the kaon-pion discrimination provided by the 
RICH system, and the absence of any neutrals in 
the final state, which allows for good reconstruc- 
tion efficiency and powerful vertex constraints. 
Consideration of the trigger and reconstruction 
efficiencies of similar topology decays reported 
in [TS] leads to the expectation of sample sizes 
of more than 1000 events per year of operation. 

3. Estimating the 7 sensitivity in 

B± ^ (K+K-7r+7r-)DK± decays 

In this section we formulate a model to describe 
B± (K+K-_7r+7r-)DK± decays. This model 
neglects D° — oscillations and CP violation in 
the D system, which is a good approximation in 
the Standard Model. The model is then used in a 
simulation study to estimate the sensitivity with 
which 7 can be determined from an analysis of 
B± (K+K-7r+7r-)DK± events. 

3.1. Decay Model 

In the same way as the kinematics of a three- 
body decay can be fully described by two vari- 
ables (Dalitz Plot), typically S12 = (pi +P2)^, 
S23 = iP2 + Ps)^, where pi,p2,P3 are the 4- 
momenta of the final state particles, so can a 
four-body decay be described by five variables. 
In this paper we use the following convention for 
labelling the particles involved in the D decay and 
their 4-momenta: 



Decay: 


D - 


^ K+ 


K- 


7r+ 




Label: 





1 


2 


3 


4 


4-mom. : 


Po 


Pi 


P2 


P3 


P4 ■ 



We also define: 



S^j EE {pi+p/f, 

Sijk = {Pi+Pj+Pkf, 
Uj EE {pi~pjf. (3) 
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We then choose a set of five variables to describe 
the decay kinematics: tgi = S234j S12, S23, S34 
and t^Q — S123. From these variables all other 
invariant masses Sij,Sijk, and, for a given frame 
of reference, all momenta pi can be calculated. 

In contrast to the phase space density for three- 
body decays, which is uniform in terms of the 
usual parameters si2,S23, four-body phase space 
density, dcj)/ dtQidsi2ds23ds3idt4Q, is not flat in 5 
dimensions, but proportional to the square-root 
of the inverse of the 4-dimensional Grahm deter- 
minant 1171: 
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(4) 



where Xij =pi-pj ^ ^ y^^j ...^ ,,„^^ 

The total decay amplitude for the decay to 
the K+K^7r+7r~ final state is the sum over all in- 
dividual amplitudes Ak to each set of intermedi- 
ate states fc, weighted by a complex factor |cfc|e*"^'' 



k 



\Ck\e 



(5) 



An analysis of the D'' K+K^tt+tt^ de- 

cay amplitude is reported in [T2] . which fits 
10 separate contributions. In this analysis, 
however, no distinction is made between the 
modes 0° Ki(1270)+K-, Ki(1400)-K+ and 
K*(892)0K-7r+, and those decays to the CP- 
conjugate final states. In our study we base \ck\ 
and (j)k on the values found in 'T2], but consider 
different scenarios for the relative contributions of 
the above modes. In order to label these scenarios 
we make the definitions 

_ |c/fi(1270) + K- P 



R 



Ki{1270)K 



and 

^</'Ki(1270) 



|c/fi(1270)-K+h 



(6) 



K = (PKi{1270) + K- 

— 4>Ki{1270)- K+ ■ (7) 

We define similar variables for the D*^ 
Ki(1400)±KT and 0° ^ K*(892)0K±7r=F decays. 



Our default scenario assumes the arbitrary values 



R 



Ki{1270)K 



= R 



Ki{1400)K 



= Rk' 



K'{892)"K7r 



= 1, 



A0_R-i(i27o)K = 39°, A(t)Ki(uoo)K = 211° and 

A4>K' {892)0 Ktt = 115°. 

The amplitudes Ak are constructed as a prod- 
uct of form factors (F;), relativistic Breit-Wigner 
functions (BW), and spin amplitudes (s/) which 
account for angular momentum conservation, 
where I is the angular momentum of the decay 
vertex. Therefore the decay amplitude with a sin- 
gle resonance is given by 



A^Fi-si-BW, 



(8) 



(where the subscript k has now been omitted), 
and a decay amplitude with two resonances a and 
P is written 



A = si-Fi^-BWo,-Fifi-BWf3. 



(9) 



For Fi we use Blatt-Weisskopf damping fac- 
tors [H] and for s; we use the Lorentz invariant 
amplitudes [H], which depend both on the spin 
of the resonance (s) and the orbital angular mo- 
mentum. 

With these definitions, the total decay ampli- 
tude for DK^,D K+K^tt+tt^ is given 

by 

A- = .4 (B- ^ (K+K-7r+7r-)j3K") 

S12, S23, S34, t4o) 



i(SB -7) 



Ad° {t02, S12, Sl4, S34, t^o) ^ . 



(10) 



The corresponding expression for the CP conju- 
gate decay is 



A+ = A(B+ ^ (K+K-7r+7r-)j^K 



Ab {^Ado ito2, Sl2: Sl4, S34j ^3o) 

(^01,Si2,S23,S34,t4o)^. (11) 
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The total probability density function for a 



(K+K- 



P- = N\A-\' 



I K event is then given by 



(12) 



dtoidsi2ds23ds3idt. 
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(with an equivalent expression for decays) 
where A'^ is an appropriate normalisation factor 
which may be obtained through numerical inte- 
gration. 

3.2. Simulation Study 

To estimate the statistical precision achievable 
with this method, we generated several simula- 
tion samples which we then fitted to determine 
the parameters of interest, most notably 7. The 
samples were generated neglecting background 
and detector effects. Figure [2] shows the pro- 
jections of the chosen kinematical variables for 
200k events, separately for B+ and B^ decays, 
and the CP-asymmetry, defined as the number of 
B"*" events minus the number of B~ events, nor- 
malised by the sum. In these projections the ob- 
servable CP violation is small, typically being at 
the few percent level only. Full sensitivity to 7 is 
obtained through a likelihood fit to all five vari- 
ables. 

A log-likelihood function is defined as: 

iog(/:) - J2 i°g(^r) + E i°g(^^)' (13) 



all B- 



all B + 



where the probability density functions are de- 
fined as in expression II 2i and the sums run over 
all B candidates in the sample. The function was 
maximised for each sample using the MINUIT 
package [I9j, with rs, Sb and 7 as the free pa- 
rameters. (It is assumed that all parameters asso- 
ciated with the D decay model are known.) Each 
sample contained 1000 events. A scan of the neg- 
ative log-likelihood, plotted for 7 against Sb is 
shown for a typical sample in Figure [3l The func- 
tion is well behaved with a minimum close to the 
input value and a second solution at 7 — 180° and 
Sb — 180°. Also shown is a scan for 7 against r^. 

For each sample the fitted parameters and the 
assigned errors were recorded. The reliability of 
the fit result was studied for each variable by 
constructing the 'pull distribution', which is the 
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Figure 2. Distributions (left) of 200k simulated 
events for the 5 kinematical variables shown for 
B+ (solid) and B" (dashed) decays separately. 
Also shown (right) are the asymmetries between 
the B+ and B^ distributions, where the asymme- 
try is defined as the number of B~ decays minus 
the number of B+ decays normalised by the sum. 
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Figure 3. Negative log-likelihood shown for 7 
against Sb (top) and 7 against rs (bottom) for 
a typical simulated experiment of 1000 events. 
The input values are at 7 = 60°, Sb = 130° and 
rs = 0.10. 



Table 1 

Result from 218 simulated experiments, showing 
the average assigned error, and the mean and 
width of the pull distribution. 



Error 



Pull mean Pull width 



7 14.4± 0.7° 
Sb 14.3± 0.5° 
rs 0.023±0.001 



0.13±0.07 
O.lOiO.07 
0.11±0.06 



1.10± 0.05 
1.09± 0.05 
0.95± 0.05 



Table 2 

Dependence of 7 fit results on the value of r^, 
showing the average assigned error and the means 
and widths of the pull distributions. 





Error 


Pull mean 


Pull width 


0.05 


24.4±0.6° 


0.12±0.21 


1.05± 0.15 


0.10 


14.4±0.7° 


0.13±0.07 


1.10± 0.05 


0.15 


8.8±0.2° 


0.00±0.29 


0.96± 0.21 


0.20 


7.2±0.1° 


-0.06±0.22 


1.11± 0.16 



difference between the fitted and input parame- 
ter, divided by the assigned error. The means 
and RMS widths of the pull distributions are dis- 
played in Table [T] and are seen to be compatible 
with and 1 respectively. This indicates that the 
log-likelihood fit is unbiased and the returned er- 
rors are reliable. The fit errors are also included 
in Table [1] averaged over all fits. 7 is extracted 
with a precision of 14°. There is very little cor- 
relation between the three fit parameters, as is 
clear from the contours in Figure [31 

The size of the interference effects in — 
DK* decays, and hence the sensitivity of the fit 
to 7, depends on the value of rs. To investigate 
this dependence several 1000 event samples were 
generated with different values of between 0.05 
and 0.20. These samples were then fitted as pre- 
viously. The fit result on 7 and associated uncer- 
tainty for each vb value are shown in Table [2] It 
can be seen that the 7 error varies approximately 
linearly with the inverse of r^. 

As explained in Section [5Tl1 the fitted model re- 
ported in [12] does not distinguish between the 
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Table 3 

Statistical uncertainty on 7 for various values of 
R and A0. These parameters are defined in ex- 
pressions [5] and [7] with the same values being used 
for Ki(1270)K, Ki(1270)K and K^(892)0K7r. 
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relative contribution of certain D decay ampli- 
tudes and their CP-conjugate final states. The 
importance of this unknown information on the 
fit sensitivity was assessed by generating and fit- 
ting 1000 event simulated datasets with different 
values of the R and parameters defined in ex- 
pressions [6] and [T] In varying these parameters 
the overall contribution of each mode and its CP- 
conjugate state, eg. |^(D° Ki(1270)+K-) + 
A(D0 ^ Ki(1270)-K+)|2, was kept constant. 
The results for the uncertainty on 7 are shown 
in Table [3] in the case where a common value of 
R and A(p is taken for the three final states under 
consideration. In a further study the phase shift 
A(j) was allowed to take different values between 
the three modes. Four scenarios were considered 
with the following arbitrary (randomly chosen) 
sets of values for A(j)K,(i27o)K, A0Xi(i4oo)k and 
A(/)K*(892)0if^ respectively: 

1. 39°, 211° and 115° (default); 

2. 53°, 108° and 15°; 



Table 4 

Statistical uncertainty on 7 for various R values 
and different A(f> scenarios. The same R values 
are being used for Ki(1270)K, Ki(1400)K and 
K*(892)"K7r. The scenarios for are given in 
the text. 



R 



A(j) scenario 



0.25 14° 19° 14° 13° 

0.50 14° 17° 14° 15° 

1.00 12° 18° 26° 20° 

2.00 13° 14° 16° 19° 



3. 55°, 344° and 173°; 

4. 209°, 339° and 87°. 

The statistical uncertainties found on 7 for these 
scenarios are given in Table [H For both Table [3] 
and Table 3] only a single experiment was per- 
formed at each point in parameter space, hence 
the stated error carries an uncertainty of a few 
degrees. However any minor variation in result 
arising from the exact value of the fitted tb pa- 
rameter, experiment-to-experiment, has been cor- 
rected for by using the dependence observed in 
the study reported in Tabled 

It can be seen that the precision of the fit 
is fairly uniform over parameter space, with a 
typical value of 15°. In certain cases however 
the precision is worse, particularly when R = 
1 and/or A(/) = 0. More detailed studies of 
D K+K~7r+7r~ decays are therefore needed 
to reliably estimate the instrinsic sensitivity of 
(K+K^7r+7r^)D for a 7 measurement. 
However, variations in other aspects of the D — > 
K+K^TT+TT^ decay structure were found to have 
limited consequences for the fit precision. 
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Finally it was investigated what biases would 
be introduced in the 7 extraction through incor- 
rect knowledge of the decay model. Experiments 
were performed in which the datasets were gener- 
ated with the full model in the default scenario, 
but fitted with a model which omitted all the de- 
cay amplitudes with a contribution less than 3% 
to the overall rate. Shifts of up to 8° were ob- 
served in the measured value of 7. This value can 
be considered as an upper bound to any final sys- 
tematic uncertainty, as it will be possible to ac- 
cumulate very large samples of D — > K+K^tt+tt^ 
events at the LHC, which will allow the decay 
model to be refined and improved with respect 
to the one assumed here. Additional information 
will also become available from CP-tagged D de- 
cays at facilities operating at the ^'(3770) reso- 
nance [20] , 

4. Conclusions 

We have shown that the decay 
(K"'"K~7r"'"7r~)D can be used to provide an in- 
teresting measurement of the unitarity triangle 
angle 7. With 1000 events and assuming a value 
of tb = 0.10 it is possible to measure 7 with a pre- 
cision of around 15° . The exact sensitivity achiev- 
able depends on the relative contributions of cer- 
tain unmeasured modes in the D decay model. 
The final state, involving only charged particles, 
and kaons in particular, is well suited to LHCb. 
A full reconstruction study is necessary to esti- 
mate reliably the expected event yields and the 
level of background. 

Finally we remark that the same technique of a 
four-body amplitude analysis in — > DK* de- 
cays can be applied to other modes, most notably 
the 'ADS' channel D ^ K^^tt^tt+tt- . 
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